The geometrically-averaged density of states as a measure of localization 



O 

o 

(N 
(N 



CD 



I 

o 
o 



> 
O 

On 
O 

O 



o 
o 



X 



Yun Song, W. A. Atkinson, R. Wortis 
Department of Physics & Astronomy, Trent University, 
1600 West Bank Dr., Peterborough ON, K9J 7B8, Canada 
(Dated: February 6, 2008) 

Motivated by current interest in disordered systems of interacting electrons, the effectiveness of tfie 
geometrically averaged density of states, Pgiio), as an order parameter for the Anderson transition 
is examined. In the context of finite-size systems we examine complications which arise from finite 
energy resolution. Furthermore we demonstrate that even in infinite systems a decline in pgiui) with 
increasing disorder strength is not uniquely associated with localization. 



I. INTRODUCTION 

Although disordered metalhc systems have been stud- 
ied extensively for many years, a reasonably complete 
theoretical understanding has only been developed for 
noninteracting particles,^ and the effect of interactions 
in disordered media is an open and challenging problem. 

In 1958, Anderson^ proposed that noninteracting elec- 
trons may be localized by the interference of wavefunc- 
tions elastically scattered off of multiple impurities even 
in situations in which the electrons are classically un- 
confined. Extensive studies^ since then have produced 
a clear picture of the Anderson transition. In d = 3 di- 
mensions, and for sufficient disorder, single-particle wave- 
functions at the band edges are localized, and decay expo- 
nentially on a scale referred to as the localization length, 
^. ^ is a function of energy, and diverges at the upper and 
lower mobility edges ujc± which separate localized and ex- 
tended states; states with energies uj > lUc+ ot lo < ojc- 
are localized, while states with lOc- < to < lUc+ are ex- 
tended. The mobility edges move towards the band cen- 
ter as disorder is increased, and a metal-insulator transi- 
tion occurs when either mobility edge crosses the Fermi 
level. 

The question of how electron-electron interactions in- 
fluence the interference effects responsible for Ander- 
son localization has received much attentionii^i^ with 
calculations using mean field, diagrammatic perturba- 
tion and renormalization techniques. While much has 
been learned, a clear understanding of the strong inter- 
action limit has not been achieved. Recently a num- 
ber of approaches have been developed to use dynam- 
ical mean field theory (DMFT) methods to study the 
Mott-Anderson metal-insulator transiton.i^'^'^'^ '^^i^^ In 
its original form DMFT is a local approximation in which 
a single site is embedded in an effective medium. The 
distinction from standard mean-field theory is that the 
on-site interaction is treated exactly, and the effective 
medium is chosen to approximate the effects of other 
interactions. One solves self-consistently for the local 
single-particle Green's function G'(r, r, z), where r is a 
position and z a complex frequency. 

DMFT captures the physics of the interaction-driven 
metal-insulator transition-"^- (the Mott transition), but 
two significant challenges arise in its application to the 



study of the Anderson transition. First, one must cap- 
ture the subtle physics of localization within an effec- 
tive medium theory and, second, one must extract infor- 
mation about nonlocal behavior from the local Green's 
function. The first problem has been the topic of recent 
work^ili^i^i^ii^ while the focus of the current article 
is the second: What quantity calculable from the local 
Green's function is a practical measure of localization? 

For numerical calculations on noninteracting-electron 
systems, a number of convenient measures exist: includ- 
ing the Thouless number^ii^ spectral statistics,-- and 
especially the inverse participation rati o^"^'^^ (IPR). A 
key point, however, is that in the absence of electron- 
electron interactions the many-body wavefunction is a 
product of well-defined single particle eigenstates, and 
all of these methods rely on knowledge of these states. 
However, in interacting systems no such states exist. 

A possible approach is suggested by the facti^iii that 
for a system with localized states the local density of 
states (LDOS), ^(r,^) = — vr^^Im G'(r, r, -I- iO+), has 
a discrete spectrum. In an infinite system the density 
of states of the whole system (DOS) is in general a con- 
tinuous function of to, and when the states are extended 
the LDOS is correspondingly continuous. However, in an 
infinite system with localized states the DOS will remain 
continuous while the LDOS becomes discretci ^^^'^^i^^ This 
has been discussed more recently in the context of non- 
interacting electrons by a number of authors fi^i^SiSii^ 

With this in mind, the geometric average of the density 
of states (GADOS) 
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has been proposed^'^'^'^i as a possible order parameter 
of the Anderson localization transition. Note that if the 
LDOS at a given energy cu is zero on any site, the GADOS 
at that energy will be zero. Therefore, in an infinite 
system Pg(0) — > signals a metal-insulator transition, 
where lu = corresponds to the Fermi level. 

In this paper we explore some potential pitfalls in the 
use of this measure. Our central results are (1) that a 
finite energy resolution changes the scaling behavior of 
Pg{uj) and makes the determination of the mobility edges 
difficult and (2) that simply a decline, as opposed to a 
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zero, in the GADOS is not uniquely associated with the 
phenomenon of locahzation. We arrive at these conclu- 
sions in the context of a system of non-interacting elec- 
trons on a lattice with random potential values at each 
site (Section III Al) . First, we consider finite-size systems 
and use a standard numerical diagonalization approach 
(Section III B[) . In these systems, the density of states is 
discrete, but a spectrum similar to the thermodynamic 
result can be obtained by inserting a finite energy reso- 
lution, a process mathematically parallel to the presence 
of an inelastic scattering rate. We demonstrate (Sec- 
tion [TlEI some misinterpretations which may arise, we 
sketch the scaling procedures, explored in more detail 
elsewhere,— which provide reliable phase information, 
and we examine in particular the influence of a finite en- 
ergy resolution on the scaling behavior. Second, we treat 
the infinite system within the coherent potential approx- 
imation (CPA) (Section III C[) . This method is known to 
accurately reproduce the disorder-averaged DOS while 
not capturing the multiple-scattering processes respon- 
sible for Anderson localization. We find (Section IIII Bp 
that the CPA GADOS declines with increasing disorder 
and suggest that this is due to small wavefunction am- 
plitudes at sites of extreme potential. In summary, the 
dependence of Pg{Lo) on disorder strength alone, in the 
absence of a clearly defined zero, is not a reliable mea- 
sure of localization, and furthermore when the energy 
resolution is finite, scaling will be affected. 



II. CALCULATIONS 
A. Model 

In this work, we will use a site-disordered model, in 
which the site energies take random values, while the 
hopping matrix element t is fixed. For the simplest case 
of spinless fermions, the Hamiltonian takes the form 

H = Y.\^eS\-tY,\^){j\ (2) 

where indicates nearest neighbour sites on the lat- 
tice, and |i) is a ket corresponding to the Wannier orbital 
attached to the ?th site. In this work, we will present re- 
sults for box-distributed disorder, where are uniformly 
distributed between —W/2 and +W/2, with probability 
distribution P(e) = W~'^Q{W/2 - |e|) where 9 is a step 
function. We will consider two cases: the d = 1 chain 
where all states are known to be localized for any nonzero 
and the d = 3 cubic lattice for which all states are 
localized when W = 16.5,^1^ 

We take two approaches to finding the electronic struc- 
ture of the disordered lattice-exact numerical calculation 
and the coherent potential approximation-both of which 
are standard. 



B. Exact numerical calculation (ENC) 

The Hamiltonian given in Eq. Q describes noninter- 
acting Fermions, so that exact single-particle eigenstates 
may be found by numerically diagonalizing, using the 
LAPACK libraries, the matrix Hij = ti5ij —t6(^ij^ , where 
S(^ijj is one for nearest neighbor sites and zero other- 
wise. For a system with N sites, we obtain eigenvalues 
Ea {a = 1, . . . N) and corresponding eigenvectors ^'^(ri). 
We consider systems of linear size L (measured in units 
of the lattice constant) and dimension d = 1 or 3, having 
N = L'^ lattice sites in total. 

To calculate the GADOS using the numerical solution, 
the LDOS is obtained from the resulting eigenvectors: 

/9(r„w) = ^|*„(r,)p5(t^-i;a). (3) 

a 

In numerical calculations, where one has a discrete spec- 
trum, one generally approximates the delta-function by 
bins of width 7: 

5{uj-E^)^^-'q(1--\E^-lo\). (4) 

The LDOS approaches the thermodynamic result when 
7 > AAr(w) where AAr(w) = [p{u)L\^ ^ Adt/N is the 
energy level spacing due to the finite size of the system. 
This binning procedure is essentially equivalent to ap- 
proximating 

p{vi,uj) w ~TT^H-mG{Yi,Yi,uj + 17), (5) 

which is significant because it provides a point of ref- 
erence between our results and results for many-body 
calculations of the single-particle Green's function. 

The DOS is the arithmetic average of this: p{lo) = 
iV^^ p(ri, w), while the GADOS is as defined in 
Eq. Representative examples of our results for the 
DOS and the GADOS in both one and three dimensions 
are shown in Fig. [TJ 

We present results both for 7 = and 7 7^ 0. For 
7 7^ 0, we calculate the LDOS using Eq. (g]). For 7 = 0, 
we calculate the LDOS using Eq. Q (a procedure which 
is only possible for noninteracting systems) in which case 
Pg is only nonzero for energies corresponding to an energy 
eigenvalue. This result is subsequently smoothed using a 
binning procedure. Note the key distinction that when 
7 = the binning is done after the GADOS is calculated. 

An important technical issue which must be discussed 
is disorder averaging. In order to improve the signal- 
to-noise ratio, we typically present results which have 
been averaged over several hundred impurity configura- 
tions. For piuj) and for the IPR, this is a meaningful 
process. For Pg{u>) it is less clear whether one should 
take the arithmetic average or the geometric average, 
as was done in Ref. [l^, of different samples. We have 
studied both averaging procedures and have determined 
that they yield similar results, except for small values of 
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FIG. 1; (Color online) (a) d = 1 DOS; (b) d = 1 GADOS; 
(c) d = 3 DOS; and (d) d = 3 GADOS. In each case results 
from the ENC and from the CPA are shown for two values 
of disorder strength. The ENC results for d = 1 are for A'' — 
2000 averaged over 400 samples, while those for d = 3 are 
for L = 15 averaged over 200 samples. All ENC results use 
7 = 0.05. 



7 Ajv, for which the geometric averaging becomes mi- 
reliable: If any particular disorder configuration fails to 
have an eigenvalue in a particular bin, then Pg{uj) van- 
ishes for that bin. For this reason, we find that geomet- 
rically averaging disorder configurations for Pg{uj) exag- 
gerates the already high sensitivity of the GADOS to 
single sites with high potential, and we have adopted the 
arithmetic average for calculating the disorder average of 

Pgi^)- 

For purposes of comparison, it is also useful to calcu- 
late the IPRfiiii^ which is known to be a reliable and 
sensitive measure of localization. For normalized eigen- 
states, the IPR for eigenstate a is 



(6) 



and it is common to average the IPR over energy bins of 
width 7 according to 
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It is the scaling of these quantities with system size which 
determines whether states are localized: For a system 
with N sites, both and A{ijj) scale as N^"^ for ex- 
tended states, but saturate at a finite value when the 
localization length ^ is less than the linear dimension L 
of the system. In practice when systems with L > S, are 
calculationally inaccessible, localization can be inferred 
when A{ll!) vs. extrapolates to a nonzero value. 



C. Coherent potential approximation (CPA) 

The CPA is the second approach we take to determine 
the electronic structure of our model. This is an effec- 
tive medium theory, modelling a randomly distributed 
impurity potential with a self-energy The self- 

consistent equation for the self-energy in the CPA is^ 
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where the on-site potential is distributed randomly ac- 
cording to the distribution P(e), for which we use 
the box-distribution as described in Section III Al The 
disorder-averaged local Green's function is 
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(9) 



where, for a tight-binding lattice 



—2t cos k 
-2i(cosfcj. -I- cos ky 



d = 1 
■ cos kz) d — 3 



is the dispersion of the lattice. Equation ((8]) is a state- 
ment that when averaged over the effective disorder po- 
tential e — S the T-matrix vanishes. 

We construct Pg{uj) from the CPA results using the 
following simple approach: We find the LDOS for a site 
embedded in an effective medium characterized by the 
CPA self-energy where is determined from 

Eq. ([5]). We call this LDOS p^{0,uj) because the site is 
defined to sit at the origin and the site potential is e. 
We then take the geometric average over an ensemble 
of systems where e satisfies the same distribution func- 
tion P(e) used in the CPA calculation: ln[p^^^(ti;)] = 
/ deP{e) ln[pe(0, w)]. Note that if we take the arithmetic 
average instead we simply recovered the CPA DOS. ie. 
/deP(e)p,(0,c^) = pCPA(^) = Gioc{u;)/n. 

Specifically, we solve the equations of motion for the 
Green's function with the site energy e at the origin, 

[a;+-e]G(0,0,c^)-t^V»>G(r„,0,c^) = 1 

71 

- I](cj)]G'(ri,rf,cj) - t^(5(i^„)G(r„,rf, w) = (5,;,^ 

n 

where the case i = £ = is excluded from the second 
equation and lu^ = uj + iO. These give 

G(0, 0, u;)^[l + Giociio)in^) - e)]-'Giociuj) 



(10) 



with Gioc{i^) defined in Eqn. We remark that if 

one imposes the self-consistency condition Giod^^) = 
/ (ieP(e)G(0, 0, w) then the procedure we have out- 
lined is equivalent to the CPA. Noting that p^{Q,ll)) = 
-Im G(0,0,tj)/7r, we then find the GADOS via 



lnp^PA(cc;) = / deP{e)\n 
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FIG. 2: (Color online) Sketches of the DOS and LDOS at two 
well-separated sites A and B in a sample of fixed size L with 
varying levels of disorder: (a) no disorder, (b) W < Wc, (c) 
W > Wc hut ^ > L, and (d) W > Wc and ^ « L and in 
particular ^ << |r,4 — rs|. 



The CPA results for the p{uj) and Pg{uj) are compared 
with exact numerical results in Fig. [TJ While the CPA 
reproduces the exact p{lu) accurately, there is a signifi- 
cant discrepancy between p'^^^luj) and Pg^'-^{uj), pre- 
sumably because the CPA does not capture the quan- 
tum interference to which the LDOS and hence GADOS 
are sensitive. Nonetheless, the CPA captures the general 
trend that Pg{u!) is a decreasing function of W. We will 
explore this point in more detail below (Section [TnB)) . 



III. RESULTS AND DISCUSSION 
A. Finite systems 

We begin with our results on finite size systems in one 
and three dimensions. It may be useful to keep in mind 
three energy scales and three corresponding length scales. 

is a measure of the average level spacing in the LDOS 
caused by localization. Larger values of correspond 
to stronger localization and shorter localization length: 
Aj ~ l/^"*. Both ^ and A^ depend on u. A^v is a 
measure of the mean level spacing in the DOS caused 
by the finite size of the system. Larger values of An 
correspond to smaller system sizes: Ajv ^ ^/N = l/L''. 
Finally, 7 is the size of the energy bin width used in 
numerical calculations of the LDOS [Eq. ^] and the 
IPR [Eq. ([7])]. As suggested by Eq. (O, 7 plays the same 
role as an inelastic scattering rate and the corresponding 
length scale is therefore an effective inelastic mean free 
path, £^ ~ 1/7!/'^,^ 

Fig. [5] presents a useful framework for considering the 
effects of finite size on calculations of pg. Four cartoons 
are presented corresponding to four values of disorder 
strength, all in a finite-size system of length L with 7 = 0: 
(a) no disorder, (b) weak disorder such that no localiza- 
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FIG. 3: (Color online) (a) The inverse participation ratio 
as a function of A'^"^ at four different energies for a d = 1 
chain with W = 1 and 7 — 0.1. (b) The localization 
length extracted by a fit to A(ijj,N) (see text), and an al- 
ternative definition {lo, Nmax) of the localization length 
{Nmax = 2000). The IPR is averaged over 400 disorder con- 
figurations for each value of A''. Energies are measured in 
units of the hopping integral t. 



tion takes place, (c) intermediate disorder such that lo- 
calization occurs but with localization length large com- 
pared to the system size, and (d) strong disorder such 
that the localization length is much less than the system 
size. In each case a sketch is given within some frequency 
window of the DOS (for the whole system) and the LDOS 
at two well-separated sites A and B. The lengths of 
the lines correspond to the weights of the delta-function 
peaks in the DOS and LDOS at the discrete eigenener- 
gies of the system. For a system with no disorder, case 
(a), the states will be not only extended but also, at low 
energies, fairly uniform so the weights of the eigenstates 
at a given site will be roughly the same. When disor- 
der is introduced but no localization occurs, case (b), 
the weights at each site are shifted as sites become more 
and less energetically favorable. When localization be- 
gins but ^ > L, case (c), the unevenness of the weights in 
the LDOS simply becomes more pronounced. Only when 
localization is strong enough that « L, case (d), do 
some of the weights become exponentially small. 

Presenting our results, we begin with the case of a 
one dimensional chain, which is useful for several rea- 
sons. First, it is simple: there is no Anderson transition. 
All states are strongly localized for any nonzero W, with 

Ri ttI where £ is the mean free pathfiiSl Furthermore, 
the limits L < ^ and L > ^ are both accessible in numer- 
ical computations. Finally, the fact that we can study 
instances of localization with W less than the bandwidth 
allows us to distinguish two distinct factors which cause 
a decline in the GADOS-exponential tails of localized 
eigenstates and small wavefunction amplitudes at sites 
of extreme potential. 

Fig. [3] demonstrates how the IPR shows localization 
in this d = 1 case. In Fig. [3] (a), there is a clear 
crossover from A ^ 1/N to A roughly constant at 
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FIG. 4: (Color online) Dependence of the GADOS at a; = 
on W in one dimension, showing ENC results for several bin 
widths and CPA results. Calculations are in the strongly- 
localized regime N — L = 2000 S> C ^^nd the results are aver- 
aged over 400 disorder configurations. Energies are measured 
in units of the hopping integral t. 



L = ^. The corresponding energy-dependent localiza- 
tion lengths are shown in Fig. [3] (b). The locahzation 
lengths have been extracted by fitting the IPR to a 
form A{uj,N) = ^(0;, c») coth(A^/2^) which we expect 
for one- dimensional wavefunctions that are decaying ex- 
ponentially over a length scale ^ with periodic bound- 
ary conditions. Because A{uj, N) saturates for L :s> 
A{uj, 00) ~ A{uj, Nmax), where Nmax = 2000. The values 
of ^ (w) determined in this way are quantitatively consis- 
tent with earlier calculations and are also close to the 
quantity A~^{ijj, 00) that is frequently taken as a defini- 
tion of the localization length. The quantitative differ- 
ences between the two definitions are not important for 
this work. Rather, the main point of Fig. [3] is to estab- 
lish that the d = 1 results presented below are for the 
strongly-locaHzed regime L :s> ^, where we expect Pg{uj) 
to be strongly influenced by the exponential tails of the 
wavef unct ions . 

Fig. |4] shows the variation of the GADOS at the band 
center as a function of disorder strength for a finite size 
system. Results from the ENC are shown for several 
values of the energy bin widths including 7 = 0, and 
the CPA result is also shown. Naively, there appears to 
be an Anderson transition at nonzero W. If it were to 
be assumed that the Anderson transition occurs when 
Pg{oj) = rip{uj) with <C 1, then Pg{u>) erroneously pre- 
dicts an Anderson transition at a nonzero critical disorder 
W^^ , where W^^ is strongly 7 dependent. This happens 
for two reasons: First, when 7 > Aj, the discrete nature 
of the LDOS is hidden, and moreover when L < ^ even 
exponentially localized states have significant amplitude 
at all sites. The resulting potential for errors points both 
to the importance of performing a proper scaling analysis 
and also to the significance of 7 in the scaling, a feature 
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FIG. 5: (Color online) (a) The scaling with chain length of 
the ENC GADOS at the band center in one dimension with 
W — 1 for several bin widths including zero. Results are 
averaged over 400 disorder configurations, (b) The same data 
as in (a), plotted agains Leff. 



key to the applicability of this method in many-body cal- 
culations. 

We remark that the 7 dependence of Pg{io) at small 
values of W in Fig. [4] is a finite size effect due to the 
varying number of discrete states which fall in each en- 
ergy bin [see Eq. ([4])] . This issue does not arise when W 
is large because in this case the eigenvalues shift signif- 
icantly from one disorder configuration to the next and 
hence the number which fall in any given energy bin is 
stabilized by sample averaging. Equivalently, finite size 
effects appear for small W because the mean free path is 
larger than the system size. 

In Fig . [5l we examine the scaling behavior of Pg(j-^)- 
Ref. [2l| has done a detailed scaling analysis of the crit- 
ical regime L < ^ for the case 7 = 6, i.e. where indi- 
vidual wavefunctions are studied. Our goal is to study 
the influence of a finite energy resolution, a necessary 
feature of a number of many-body approaches, on the 
scaling behavior. Figure [5^ a) shows the scaling of /9g(0) 
with chain length for Vl^ = 1 and several values of 7 
including zero. At this disorder strength, ^(0) « 120 
lattice spacings and Aj(0) « 0.05i. For 7 = and 
L > ^, as the system size increases the exponentially 
decaying form of the localized states manifests itself, and 
Pg(0) p(0) exp(— L/2^). (This form is derived using 
periodic boundary conditions). For 7 = 6 and L < ^, 
the system is in a critical regime where the decline has a 
power law rather than exponential form>Si however the 
difference is difficult to distinguish over this short scale. 
For nonzero 7, pg scales exponentially for L < £j but 
saturates at larger values of L. We make the ansatz that 
there is an effective system size dominated by the shorter 
of the two relevant lengths: 

L,ff = {l/L + l/l,r\ (12) 
where l-y = (pi)^^'^, such that Pg(0) ~ 
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FIG. 6: (Color online) Dependence of the GADOS at the 
band center on W (in units of t) in three dimensions showing 
ENC results for several bin widths and CPA results. Inset 
includes results for weak disorder where there is a strong 7 
dependence. L = 15, A'' = L"^ and the results are averaged 
over 200 disorder configurations. 



p(0) exp{—Lgff /2^). In Fig. ^h), a plot of the data from 
Fig. inja) versus L^g shows that the ansatz works well. 
One immediate consequence of Eq. (jl2p is that when 



L :$> £-y, Pg{0) saturates at ~ p(0) exp(— £^/2^). Thus, 
in Fig. m the apparent critical disorder W^^ actually 
marks the crossover from £^ > £j to < £j rather than 
a localization transition. We conclude that the scaling 
of Pgiio) with N can be used to place an upper bound 
on the disorder strength at which localization occurs: If 
an exponential decay of Pg(0) with increasing N is seen, 
then localization is occurring. However, if 7 is too large, 
the effect will be masked. Essentially this analysis may 
lead to false negatives but not to false positives, and in 
particular the true result is approached as 7 — ^ 0. 

We turn our attention now to the d — i cube with 
N = sites, and L < 20. This case differs from d = 1 
in two important ways. First, there is known to be an 
Anderson transition which, for uj = occurs at Wc = 
16. 5t. Second, unless W is very large, we are generally 
restricted to L ^ ^. 

Consider first the variation of the GADOS at the band 
center as a function of disorder strength. Fig. [51 As in the 
d = 1 case, the ENC results are 7 dependent, increasing 
with increasing 7. Unlike in the d = I case, they are 
similar to the CPA result, a point discussed further in 
Section IIII Bl The weak disorder results are included in 
the inset and show strong 7 dependence for the same 
reasons discussed in the d — 1 case. 

Again, the correct approach is a scahng analysis. Fig. [7] 
(a) shows the logarithm of the GADOS versus the sys- 
tem size L {N = L^) using 7 = 0, and Fig. [7] (b) shows 




-2 



-10 



-12 



(b) -pO.OS 



■pO.O 
■pO.05 



4 6 8 



10 



15 



20 



FIG. 7: (Color online) The scaling with system size of the 
ENC GADOS at the band center in three dimensions, showing 
disorder strengths above and below the localization transition, 
Wc = 16.5. (a) 7 = and (b) 7 = 0.05. The inset shows the 
scaling with L^g for W = 20. Results are averaged over 200 
disorder configurations. Energies are measured in units of the 
hopping integral t. 



the same quantity for a nonzero energy bin width. From 
these plots it is easy to distinguish localized from dclo- 
calized behavior, at least outside of the critical regime 
L < In the delocalized phase, there is no signifi- 
cant variation with system size. This is consistent with 
the log-normal distribution of the LDOS predicted by 
Mirlinj^Sii^ and furthermore the decline in the value of 
Pg(0) with increasing W is consistent with a decrease in 
the position of the maximum of the distribution as the 
transition is approached, also predicted by Mirlin.^"-^'^ In 
the localized phase, the exponential decay as well as the 
saturation for 7^0 are similar to what is seen in one 
dimension. The larger error bars on the 7 = 0, = 16 
results are related to the proximity to the transition: The 
mobility edges are very close to w = and hence the 
energy window used in generating Pg^°(0) must be espe- 
cially narrow to avoid including localized states outside 
the mobility edge. The reduction in the number of states 
included leads to noisier results. 

In one dimension, we showed that the strong depen- 
dence of Pg{(^) on 7 can be understood in terms of an 
effective length L^jj- We find this to be true for d = 3 as 
well. The inset to Fig. ^h) shows that for the strongly- 
localized phase {W = 20), the data for 7 = and 
7 = 0.05 collapse onto a single curve when plotted against 
Leff [Eq. dm)]. We remark that for 7 = 0.05, £^ « 8 lat- 
tice constants. Since L^jj < £-y, this means that even for 
a large range of L, the effective range sampled by scaling 
is quite small. In particular, if this re-scaling of lengths 
applies to the critical regime, this implies that the deter- 
mination of the mobility edges by scaling analysis may be 
severly limited in, for example, many-body calculations 
which are performed at complex frequencies to + ij. 

Finally, Fig. [5] compares three phase diagrams for 
d — 3. First, we have shown published?! results for the 
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FIG. 8: (Color online) The bandwidth Wband and the mobility 
edges ijJc± in three dimensions, with the latter estimated using 
the scaling of the IPR and Pg{io) < ??p(a;) for two values of 77. 
Further details provided in the text. 

mobility edges based on the scaling of the IPR as a stan- 
dard of comparison. The second method used is to call 
states at a given energy localized when the curve of Pg{uj) 
versus W, like that in Fig. [51 drops below some cutoff 
Pg{uj) = r]p{uj). (System size L = 15 was used, averaged 
over 100 to 200 samples, with 7 — 0.05.) We show curves 
using two values of 77: one which results in a localization 
transition at smaller W than the IPR transition, and an- 
other which results in a localization transition at higher 
disorder strength than the IPR transition. The scaling 
of Pg{uj) with system size could also be used to produce 
a phase diagram, and the transition would always lie on 
the high W side of the IPR line, with the two coinciding 
in the limit 7 — > 0. However, the construction of such 
a diagram is similarly computationally intensive to the 
equivalent IPR calculation and in an interacting system 
offers no obvious advantages. 

B. Infinite systems 

We now turn our attention to infinite systems. Re- 
cent DMFT worls^iii^ includes plots similar to those in 



Fig.[6l showing the GADOS declining as a function of W . 
Here, we point out that factors other than exponentially 
decaying tails of localized states may lead to a decline 
(although not a zero) in the GADOS. 

The GADOS is very sensitive to small values of the 
LDOS on individual sites. Small values may occur due 
to localization, but they may also occur due to the re- 
duction in amplitude of extended wavefunctions at sites 
where the magnitude of the potential is large. The CPA 
is a convenient tool for examining this issue because it 
is known to include the latter effect but not the former. 
We show the CPA results for Pg{0) as a function of W 
for d = 1 in Fig. [5] and for c? = 3 in Fig. [51 In both 
cases p^^^(O) dechnes at large W . In the d — \ case, 
the fact that Pg^^(O) is roughly constant over the range 
< W < 4 suggests that the suppression of pg in the 
ENC calculations is due primarily to the exponentially 
decaying tails of the localized eigenstates. For c? = 3 
the decline in p^^^(O) is especially steep and fairly close 
to the ENC results, suggesting that in contrast to the 
d — 1 case the W dependence of Pg{uj) is dominated by 
lattice sites with extreme potential, rather than by expo- 
nential tails of localized eigenstates. Therefore, even in 
an infinite system, a decline in Pg{uj) with increasing W 
does not necessarily imply the approach to a localization 
transition. This does not appear to change the central 
conclusions of Refs. 001111, which are based on zeros in 
the GADOS. 



IV. CONCLUSIONS 

In summary, for finite-size systems we have demon- 
strated some of the false conclusions which may be drawn 
if a scaling analysis is not used and shown how the use of 
a finite energy resolution will influence the scaling. More- 
over, in infinite systems we have demonstrated that the 
GADOS can decline as a function of disorder strength 
without localization occurring. 
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